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Abstract 


The  results  of  research  sponsored  by  this  grant  address  two  fundamental  questions  relevant  to  the 
advancement  of  THz-regime  device  technology  (0.1  -  3  THz).  First,  they  establish  a  detailed, 
comprehensive  fundamental  understanding,  of  THz-regime  conductivity,  including  the  effects  of 
materials  imperfections  (impurities,  defects,  surface  roughness)  in  conductive  materials 
(including  metals,  semiconductors,  and  graphene).  This  includes  providing  accurate  predictive 
theoretical  models,  robustly  validated  by  experimental  measurements.  Second,  they  establish  a 
fundamental  understanding  of  low-emission-barrier  scandate  cathodes  and  identify  related, 
alternative  cathode  materials  systems  for  advanced  vacuum  electronic  cathodes  for  high  power 
THz-regime  sources. 


Summary  of  research  results 


The  original  objectives  of  this  research  grant  were: 

•  to  establish  a  fundamental  understanding  of  the  terahertz-regime  (0.1  -  3  THz)  effective 
conductivity  of  conductive  materials,  incorporating  realistic  non- idealities,  and 

•  to  establish  a  fundamental  understanding  of  low-emission-barrier  scandate  cathodes  and 
explore  related,  alternative  cathode  materials  systems. 

Both  objectives  were  met  and  exceeded  beyond  our  original  expectations.  Key  accomplishments 
included: 

1 .  We  successfully  fabricated  high  purity  copper  film  samples  with  large  (~  1  pm)  grains  and 
controlled  (varied)  surface  roughness  (15-75  nm).  Using  these  samples  we  solved  the 
challenge  of  determining  a  quantitatively-accurate  predictive  model  for  THz-regime 
conductivity  including  the  effects  of  surface  roughness.  The  details  of  this  work  have  been 
submitted  for  publication  [3]  and  are  also  published  in  Matt  Kirley’s  Ph.D.  dissertation, 
pending  approval  by  his  defense  committee  in  December  2014. 

a.  With  high  purity  (low  impurity  atom  concentration)  and  grain  sizes  much  larger 
than  electron  scattering  mean  free  path,  we  were  able  to  isolate  the  effect  of 
surface  roughness  on  the  THz-regime  conductivity  of  copper,  separating  surface- 
roughness  scattering  from  impurity  or  grain-boundary  scattering  on  ohmic 
dissipation. 

b.  We  successfully  measured  the  DC  and  THz-regime  conductivity  of  copper  film 
samples  at  400,  650,  and  850  GHz  with  varied  surface  roughness. 

c.  We  successfully  compared  the  measured  results  with  predictions  of  numerous 
theoretical  models  and  were  able  to  conclude  that  combining  measurements  of 
DC  conductivity,  AFM  measurements  of  surface  roughness,  and  application  of  the 
Hammerstad-Bekkadal  model  for  high-frequency  conductivity  including  surface- 
roughness  [1]  provides  a  reliable,  quantitatively-accurate  predictive  model  for 
THz-regime  conductivity. 


Fig.  1.  Projection  of  AFM  measurement  of  a  copper  sample  after  roughening  by  acid  etch. 


Fig.  2.  THz  conductivity  measurements  and  predictions  for  copper  surfaces  with  rms  roughness  height  equal  to  a) 
16.1  nm,  b)  60.1  nm,  and  c)  73.1  nm.  The  data  confirm  that  the  Hammerstad-Bekkedal  (HB)  model  is  a  reliable 
predictor  of  the  conductivity  including  the  effects  of  surface  roughness.  They  also  show  that  a  more  sophisticated, 
second-order -perturbation-theory  (SPM2)  model  is  a  less  accurate  predictor  of  the  conductivity.  This  means  that 
future  design  of  THz-regime  metallized  components  can  be  done  reliably,  without  need  for  guestimating  the  surface 
conductivity  of  the  surfaces. 


IQ  ^ 

2.  We  successfully  fabricated  highly-doped  (n  =  6  x  10  cm' )  silicon  samples  and 
measured  their  THz-regime  conductivity  (0.4  -  0.85  GHz).  This  research  demonstrated 
and  explained  the  difference  between  THz-regime  conductivity  of  moderately-doped  Si 
and  highly-doped  Si.  We  are  drafting  a  manuscript  of  this  work  for  publication  and  it  is 
also  summarized  in  an  appendix  of  Matt  Kirley’s  PhD  dissertation,  pending  approval  by 
his  defense  committee  in  December  2014. 

a.  Our  prior  measurements  of  moderately-doped  Si  ((n  =  1012-1014  cm'3)  (supported 
with  a  prior  grant)  demonstrated  that  conventional  Drude  theory  does  not  suffice 
to  accurately  predict  THz-regime  conductivity  [2],  In  contrast,  under  this  grant  we 


have  determined  that  conventional  Drude  theory  does  suffice  as  a  reliable, 
quantitative  predictor  of  THz  conductivity  for  highly-doped  Si. 

b.  We  have  been  able  to  explain  these  discrepancies.  At  moderate  doping  densities 
(n  <  1017),  the  quasi-static  assumption  of  conventional  Drude  theory,  i.e.,  that  cot 
«  1,  does  not  hold.  Furthermore,  at  moderate  and  low  doping  densities,  carriers 
are  scattered  mainly  by  electron-phonon  interactions  and  the  conventional  Drude 
theory’s  assumption  of  a  single  relaxation  time,  r,  is  not  representative  of  the 
distribution  of  scattering  times  experienced  by  charge  carriers.  At  high  doping 
density,  however,  ( n  >  1017  cm'3),  the  scattering  is  much  more  frequent,  such  that 
cot  «  1.  Furthermore,  at  high  doping  density,  all  interactions  become 
increasingly  electrostatically  screened  while  a  greater  number  of  ionized 
impurities  are  present  for  scattering  electrons.  Thus,  scattering,  dominated  by 
Coulomb  collisions,  is  weakly  dependent  on  electron  energy  so  that  a  single  value 
of  t  is  a  better  approximation  to  the  true  r  distribution  than  in  the  moderate  and 
low-doping  case. 

3.  We  developed  a  self-consistent  multiphysics  tool  for  examining  the  carrier  dynamics 
of  two-dimensional  electron  systems  and  applied  it  to  calculating  the  THz 
conductivity  of  graphene  supported  on  a  substrate.  The  details  of  this  computational 
tool  were  published  in  [4]  and  part  of  Chapter  3  in  Nishant  Sule’s  PhD  dissertation. 

a.  The  multiphysics  computational  tool  comprises  ensemble  Monte  Carlo  (EMC),  a 
semi-classical  technique  used  to  numerically  solve  the  Boltzmann  transport 
equation  in  order  to  simulate  carrier  motion,  coupled  with  finite-difference  time- 
domain  (FDTD)  and  molecular  dynamics  (MD).  FDTD  method,  used  for  solving 
the  Maxwell’s  curl  equations  and  MD,  used  for  calculating  the  inter-particle 
electrostatic  force,  are  used  for  coupling  electrodynamics  with  the  carrier  motion 
in  a  self-consistent  manner. 

b.  The  EMC/FDTD/MD  method  is  extended  for  simulating  carrier  dynamics  in 
supported  graphene  systems  with  charged  impurities  present  near  the  interface  of 
graphene  and  the  substrate. 


Air 


Impurities 

Fig.  3.  Schematic  of  the  simulated  structure:  single-layer  graphene  rests  on  an  Si02  substrate,  with  charged 
impurities  present  near  the  interface  between  the  two.  Carrier  transport  is  simulated  by  including  both  electrons  and 
holes  in  the  graphene  layer,  while  the  positively  charged  ions  near  the  interface  and  within  the  Si02  substrate  remain 
stationary. 


4.  We  examined  the  effect  of  charged  impurity  distributions  on  carrier  transport  in  graphene 
on  SiC>2  using  the  EMC/FDTD/MD  method.  The  results  of  this  work  were  published  in 
[5]  and  are  contained  in  Chapter  4  of  the  PhD  dissertation  of  N.  Sule. 

a.  We  showed  that  clustered  impurities  of  characteristic  width  40-50  nm 
generate  20-nm-sized  electron-hole  puddles  similar  to  experimental 
observations.  Moreover,  clustered  impurities  play  an  important  role  in  the 
carrier  dynamics  of  supported  graphene. 

b.  The  carrier  density  dependent  THz  conductivity  o(n)  calculated  using 
clustered  impurity  distributions  reproduced  the  important  features  observed 
in  experimental  measurements.  For  example:  The  residual  conductivity  and  the 
linear-region  slope  of  the  conductivity  versus  carrier  density  dependence  are 
determined  by  the  impurity  distribution  and  the  measured  slope  can  be  used  to 
estimate  the  impurity  density  in  experiment.  Furthermore,  the  high-density 
sublinearity  in  the  conductivity  stems  from  carrier-carrier  interactions. 
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Fig.  4.  Carrier  density  distribution  (blue,  electrons;  red,  holes)  depicting  the  electron-hole  puddles  formed 
in  graphene  at  the  Dirac  point  for  (a)  uniform  random  (kc  =  22  nm)  and  (b)  clustered  (kc  =  46  nm)  impurity 
distributions,  both  with  impurity  sheet  density  equal  to  5xlOn  cm”2.  The  average  size  of  the  electron-hole 
puddles,  kp,  is  estimated  from  the  FWHM  (yellow  ring)  of  the  normalized  carrier  density  spatial 
autocorrelation  function,  shown  in  (c)  and  (d),  corresponding  to  the  random  and  clustered  impurity 
distributions  from  (a)  and  (b),  respectively.  The  estimated  puddle  size  from  (c)  is  kp  =  6  nm  and  that  from 
(d)  is  kp  =  20  nm.  (e)  Characteristic  electron-hole  puddle  size  kp  as  a  function  of  the  average  impurity 
cluster  size  kc.  Each  data  point  corresponds  to  a  single  value  of  the  clustering  parameter  Lc  (swept  from  0  to 
60  nm  in  the  increments  of  5  nm)  and  is  the  average  of  14  simulation  runs;  the  error  bars  denote  the 
standard  deviations.  Insets:  Illustrative  impurity  distributions,  nearly  random  on  the  left  (Lc  =  10  nm)  and 
highly  clustered  on  the  right  (Lc  =  50  nm). 
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Fig.  5.  Conductivity  of  graphene  on  Si02  for  (a)  uniform  random  ( Lc  =  0,  Xc  =  22  nm)  and  (b)  clustered  (Lc  =  50  nm, 
Xc  =  46  nm)  impurity  distributions,  at  impurity  sheet  densities  of  1011  cm”2  (triangles),  5xl0n  cm”2  (squares),  1012 
cm”2  (diamonds),  and  without  impurities  (circles). 
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Fig.  6.  Inverse  slope  of  o(n)  as  a  function  of  the  sheet  impurity  density  for  graphene  on  Si02  at  room  temperature. 
Squares  denote  the  uniform  random  impurity  distribution  (Lc  =  0,  Xc  =  22  nm),  while  triangles  correspond  to 
clustered  impurity  distributions  that  would  give  realistic  electron-hole  puddle  sizes  ( Lc  =  50  nm,  Xc  =  46  nm).  The 
horizontal  lines  a-d  correspond  to  the  inverse  slope  values  obtained  in  several  experiments:  line  a,  L.  Vicarelli  et  al., 
Nat.  Mater.  11,  865  (2012);  line  b,  F.  Schedin  et  al.,  Nat.  Mater.  6,  652  (2007);  line  c,  Novoselov  et  al.,  PNAS 102, 
10451  (2005);  and  line  d,  Novoselov  et  al.,  Science  306,  666  (2004).  The  impurity  density  range  between  the 
intercepts  of  an  inverse-slope  horizontal  line  with  the  clustered  and  random  distribution  curves  (i.e.,  the  range  within 
the  lightly  shaded  area)  yields  an  estimate  of  the  impurity  density  range. 


5.  We  examined  the  room-temperature  complex  conductivity  o(co)  of  suspended  and 
supported  graphene  at  terahertz  frequencies  (100  GHz-10  THz)  for  a  wide  range  of 
electron  (n  =  1012— 1013  cm"2)  and  impurity  densities  (Ni  =  8xl010-2xl012  cm"2)  by 


employing  the  EMC/FDTD/MD  method.  The  results  of  this  work  were  published  in  [6] 
and  are  contained  in  Chapter  5  of  the  PhD  dissertation  ofNishant  Sule. 

a.  We  obtain  excellent  agreement  between  our  calculation  with  clustered 
impurities  and  the  experimentally  measured  o(co),  while  the  o(co)  for  a  uniform 
random  distribution  of  impurities  does  not  yield  good  agreement  for  any 
reasonable  impurity  density.  Thus,  providing  more  evidence  to  the  important 
role  of  clustered  impurity  distributions  on  carrier  transport  in  supported 
graphene. 

b.  We  show  that  the  choice  of  the  substrate  (Si02  or  h-BN),  which  is  relevant  for 
modifying  the  conductivity  with  substrate-engineering  is  important  only  at 
frequencies  below  4  THz  and  when  N;/n  <0.1.  Electron-impurity  interactions 
dominate  for  N;/n  >0.1,  and  transport  enters  the  electron-hole  puddle  regime  for 
N;/n  >  0.5. 

c.  We  find  that  the  simple  Drude  model,  with  an  effective  scattering  rate  r  and 
Drude  weight  D  as  parameters,  fits  the  calculated  <t(oj)  for  supported  graphene 
very  well,  owing  to  electron-impurity  scattering,  r  decreases  with  increasing  n 

_ I/O 

faster  than  n  and  is  insensitive  to  electron-electron  interaction.  Both  electron- 
electron  and  electron-impurity  interactions  reduce  the  Drude  weight  D ,  and  its 
dependence  on  n  is  sublinear. 


Fig.  7..  Real  part  of  g(co)  as  a  function  of  frequency.  Calculated  values  based  on  the  EMC/FDTD/MD 
simulation  (squares,  triangles,  and  circles)  match,  respectively,  the  experimental  data  from  N.  Rouhi  et  al., 
Nano  Res.  5,  667  (2012)  (blue:  UC  Irvine),  L.  Ju  et  al.,  Nat.  Nanotechnol.  6,  630  (2012)  (red:  UC 
Berkeley),  and  L.  Ren  et  al.,  Nano  Lett.  12,  3711  (2012)  (black:  Rice  University).  Dashed  lines  are  Drude- 
model  fits  to  the  simulation  data. 


Fig.  8.  Generalized-Drude-model  best  fits  to  the  calculated  a(o>)  for  graphene  on  Si02  with  different 
impurity  densities,  from  top  (red)  to  bottom  (purple):  N;  =  0,  8X1010,  1011,  2x10' ',  4xl0n,  6x10n,  8x10n  , 
1012 ,  and  2x  1012  cm-2 .  Inset:  T  corresponding  to  the  generalized-Drude-model  fits  in  the  main  panel  as  a 
function  of  N;.  Circle  colors  in  the  inset  correspond  to  the  curve  colors  in  the  main  panel. 


6.  Using  Density  Functional  Theory  (DFT)  we  have  successfully  explained  the 
mechanism  by  which  scandium  oxide  SC2O3,  (or  scandate)  cathodes  achieve  superior 
cathode  emitter  properties.  While  alternative  theories  have  been  ventured,  our 
explanation  is  the  simplest,  is  consistent  with  familiar  surface  dipole  and  impurity  doping 
mechanisms  for  lowering  work  functions,  and  explains  how  SC2O3  conducts  current  at 
high  temperatures  sufficiently  well  to  function  as  a  copious  thermionic  emitter  of 
electrons  with  good  long-life  stability  (due  to  low  volatility  of  the  SC2O3)  .  These  results 
have  been  published  in  [7]  and  [8]. 

a.  The  principal  results  of  [7]  include 

i.  Density  Functional  Theory  (DFT)  used  to  investigate  intrinsic  defect 
formation  as  a  mechanism  by  which  bulk  SC2O3  (an  insulator)  could  act  as 
a  good  conductor  in  the  thermionic  cathode  environment  (T  ~  1200  K,  P  ~ 
10'10  Torr)  to  obtain  a  better  understanding  of  observed  enhanced  emission 
from  scandate  thermionic  devices. 

ii.  From  calculating  the  formation  energies  of  relevant  intrinsic  defects  in  the 
bulk  SC2O3  structure,  and  further  by  developing  a  defect  model  whereby 
the  defect,  electron  and  hole  concentrations  were  calculated  by 
maintaining  a  self-consistent  value  of  the  Fermi  energy,  it  is  expected  that 
SC2O3  will  behave  natively  as  a  p-type  semiconductor  under  cathode 
operating  conditions. 

iii.  By  using  a  simple  Ohm’s  law  argument,  it  was  found  that  to  maintain  the 
experimentally  observed  emission  current  densities  of  approximately  10 
A/cm  through  a  200  nm  SC2O3  film,  the  SC2O3  must  have  a  conductivity 
of  1 0‘4  O'1  cm'1  or  higher  to  avoid  dielectric  breakdown.  A  very  small 


impurity  concentration  (order  of  1  part  per  billion)  is  necessary  to  realize 
this  minimum  conductivity  value  because  intrinsic  defects  alone  do  not 
produce  a  sufficient  number  of  free  carriers,  and  these  impurities  are 
expected  to  be  present  in  all  experimental  SC2O3  samples. 
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Fig.  9.  Sc203  cubic  bixbyite  conventional  unit  cell  depicting  native  and  interstitial  Wyckoff  positions  used 
in  this  study. 


Fig.  10.  Calculated  defect  formation  energies  for  intrinsic  defects  in  Sc203  as  a  function  of  Fermi  level  for 
(left)  GGA  and  (right)  HSE  functionals.  Intermediate  values  of  jusscc2°3  and  juSq2°3  (cathode  operating 
conditions:  T=  1200  K  and  P  =  10'10  Torr)  were  used  in  both  cases.  The  shaded  portion  indicates  Fermi 
energies  that  lie  above  the  calculated  GGA  bandgap. 


b.  The  principal  results  of  [8]  include: 

i.  DFT  was  used  to  calculate  the  surface  energy  barrier  (work  function)  and 
thermodynamic  stability  of  varying  coverages  of  atomic  Ba  and  Ba-0 
dimer  adsorption  on  the  (Oil)  and  (111)  SC2O3  surface  terminations  to 
investigate  the  role  Ba  plays  in  lowering  the  work  function  of  SC2O3  in 
thermionic  cathodes. 

ii.  It  was  found  that  for  bare  SC2O3  surface  terminations,  the  order  of  stability 
(from  most  to  least  stable)  is:  (111)  >  (Oil)  >  (001).  Additionally,  it  was 
found  that  atomic  Ba  adsorption  was  unstable  for  all  surface  coverages 
and  terminations  considered.  However,  partial  monolayer  coverages  of 
Ba-0  on  SC2O3  (Oil)  were  found  to  be  stable  relative  to  the  chosen 
reference  state  of  BaO.  These  stability  calculations  indicate  that  real 
devices  containing  SC2O3  will  contain  mainly  (111)  and  (01  l)-oriented 
surface  facets,  and  Ba  will  preferentially  adsorb  on  the  (Oil)  terminations 
as  a  sub-monolayer  of  BaO. 

iii.  The  lowest  calculated  surface  energy  barriers  for  electron  removal  when 
atomic  Ba  is  present  were  2.12  and  2.04  eV  for  the  (Oil)  and  (111) 
surfaces,  respectively.  The  lowest  surface  energy  barrier  for  Ba-0 
coverage  was  1.21  eV  on  SC2O3  (01 1)  for  a  partial  monolayer.  This  surface 
configuration  provides  a  plausible  and  straightforward  explanation  of 
observed  enhanced  emission  in  scandate  thermionic  devices. 

iv.  The  mechanism  of  surface  barrier  lowering  was  found  to  consist  of  a 
surface  dipole  component  as  well  as  an  electron  (n-type)  doping 
component.  It  was  found  that  combined  n-type  doping  and  surface  dipole 
modification  yielding  a  surface  barrier  lowering  that  was  less  than  the  two 
contributions  would  yield  separately. 


Fig.  11.  Sc203  (Oil)  and  (111)  slabs  depicting  fully  relaxed  Ba-O  dimer  geometries  for  several  example 
coverages.  (A),  (B)  and  (C)  are  coverages  of  1,  3,  and  7  Ba-O  per  (011)  surface  unit  cell.  (D),  (E)  and  (F) 
are  coverages  of  1,  6  and  9  Ba-O  per  (111)  surface  unit  cell.  A  stable  chain  geometry  of  O-Ba-O  bonds 
forms  at  coverages  between  3-6  Ba-O. 


Fig.  12.  Surface  barrier  as  a  function  of  Ba-0  coverage  for  the  (Oil)  termination  (A)  and  (111) 
termination  (B).  Solid  blue  lines  are  surface  barriers  of  the  Ba-O-terminated  surface  while  dashed  red 
lines  are  the  surface  barriers  of  the  opposing  bare  surfaces.  Surface  energy  as  a  function  of  Ba-0  coverage 
for  the  (Oil)  termination  (C)  and  (111)  termination  (D)  are  plotted.  The  dashed  lines  in  (C)  and  (D)  mark 
the  surface  energy  of  the  bare  reference  slabs  for  comparison.  Stable  dimer  coverages  (relative  to  BaO)  on 
(Oil)  are  circled,  characterized  by  surface  energies  which  are  lower  than  the  bare  reference  (Oil)  surface 
energy. 


7.  Using  Density  Functional  Theory  (DFT)  we  have  successfully  identified  a  brand-new 
class  of  materials  based  upon  the  perovskite  structure  that  have  potential  to  provide 
superior  high  power  microwave  (vacuum  electronic)  device  cathodes  (thermionic  or  field 
emission).  This  is  the  first  instance  of  identifying  an  advanced,  superior  cathode 
material  candidate  using  ab  initio ,  quantitative  computational  modeling.  Future 
research  using  high-throughput  screening  of  various  perovskite  alloys  has  the  potential  to 
identify  candidate  compounds  that  would  provide  lower  work  function,  longer  lifetime 
cathodes  without  the  need  for  volatile  impregnates  such  as  barium.  Our  results  under  this 
grant  have  made  such  high-throughput  computational  screening  possible  by  identifying  a 
bulk  alloy  property  that  is  an  approximate  reliable  predictor  of  low  work  function.  This 
bulk  property  is  much  faster  to  calculate  than  surface  properties  such  as  work  function. 
There  are  many  important  benefits  that  could  be  realized  from  such  work,  including 
low-cost,  long-lifetime  HPM  devices,  and  longer  lifetime  for  critical  defense 
technologies  such  as  high  power  radar  and  ECM  transmitters  and  vacuum 
electronic  (traveling  wave  tube)  satellite  power  amplifiers.  A  manuscript  of  these 
results  is  being  drafted  for  publication.  The  principal  results  include: 

a.  DFT  was  used  to  calculate  the  (OOl)-oriented  AO  and  B02-terminated  work 
functions  of  eighteen  different  technologically  relevant  3d  transition  metal 
perovskites:  the  LaB03  series  (B  =  Sc,  Ti,  Y,  Cr,  Mn,  Fe,  Co,  Ni),  SrB03  series 
(B  =  Ti,  V,  Fe,  Co),  Lai.xSrxMn03  (x=  0.0625,  0.125,  0.25,  and  0.375)  (LSM), 


LaA103  and  Bao.5Sro.5Coo.75Feo.25O3  (BSCF).  The  purpose  of  investigating  these 
materials  was  threefold:  to  ascertain  the  range  of  work  function  values  these 
materials  can  exhibit;  to  understand  the  physics  of  what  governs  the  work 
function  of  perovskite  materials,  and  to  find  a  new,  low  work  function  material 
for  application  in  high  power  devices  and  thermionic  energy  conversion. 

b.  The  work  function  range  of  these  materials  is  broad  and  varies  from  as  low  as 
1.79  eV  for  AO-terminated  SrV03  to  6.87  eV  for  B02-terminated  LaA103.  SrV03 
is  the  lowest  work  function  material  examined  here,  and  presents  itself  as  a 
promising  candidate  to  replace  existing  thermionic  cathode  materials  given  its  low 
work  function,  metallic  conductivity,  and  stable  (001)  surface  terminations. 
Additionally,  it  was  found  that  doping  SrV03  with  Ba  may  result  in  a  Ba-rich 
surface  that  could  exhibit  a  work  function  as  low  as  1.07  eV. 

c.  We  utilized  the  O  2 p  band  center  calculated  from  the  electronic  density  of  states 
of  the  bulk  materials  as  an  electronic  structure  descriptor  to  develop  an 
understanding  of  what  governs  the  work  function  physics,  and  found  there  is 
semi-quantitative  linear  correlation  between  both  the  AO  and  BO2  work  functions 
and  the  value  of  the  O  2 p  band  center  of  the  corresponding  bulk  materials.  This 
relationship  may  be  used  in  future  work  of  high-throughput  computational 
screening  of  the  O  2 p  band  centers  of  many  bulk  perovskite  materials  to  survey  a 
larger  composition  space  for  stable  materials  that  are  expected  to  fall  within  an 
expected  work  function  range,  e.g.  searching  for  other  low  work  function 
materials. 
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Fig.  13.  Trend  of  (001)  AO  and  B02  terminated  surface  work  functions  for  the  18  perovskite 
materials  as  a  function  of  B-site  element  across  the  periodic  table.  The  solid  (open)  symbols 
connected  with  a  solid  (dashed)  line  are  the  B02  (AO)  work  functions,  respectively.  Red,  blue, 
purple,  and  green  symbols  signify  the  LaB03,  SrB03,  LSM,  and  BSCF  series  of  materials, 
respectively.  LaA103  is  not  shown  on  the  plot  but  has  calculated  work  functions  of  3.25  eV  and 
6.87  eV  for  the  AO  and  B02  terminated  surfaces,  respectively. 


Fig.  14.  Plots  of  calculated  work  functions  for  the  AO  (A)  and  B02-terminated  surfaces  (B)  of 
AB03  materials  as  a  function  of  the  O  2 p  band  center  of  bulk  ABO3  materials.  In  both  plots,  the 
blue  symbols  represent  insulating  perovskites  while  the  red  symbols  represent  metallic 
perovskites.  In  (A)  there  is  a  semi-quantitative  linear  correlation  for  AO  work  function  versus 
bulk  O  2 p  band  center.  In  (B)  there  is  also  a  linear  correlation  of  B02  work  function  with  the  bulk 
O  2 p  band  center. 
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